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Excitonic versus electron-hole liquid phases in Tm[Se,Te] compounds 
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We discuss, from a theoretical point of view, excitonic phases at the pressure induced 
semiconductor-semimetal transition in TmSeo.45Teo.55, focusing, in particular, on the stability 
against an electron-hole liquid. The electron-hole pair density parameter rg[Eg, T) is calculated 
within the quasi-static plasmon pole approximation as a function of temperature T and energy 
gap Eg and converted into —Eg{rs,T). A comparison of this quantity, which is the electron- 
hole pair chemical potential, with the exciton binding energy reveals that excitons should be 
suppressed in contrast to experimental evidence for excitonic phases. We suspect therefore 
inter-valley exciton scattering and exciton-phonon scattering to substantially stabilise excitons 
in TmSeo.45Teo.55- 
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1. Introduction 

The possibility of an excitonic insulator - understood 
as a macroscopic, phase coherent quantum state (exciton 
condensate) - separating below a critical temperature a 
semiconductor (SC) from a semimetal (SM), has been 
anticipated by theorists long time ago.^ Experimental ef- 
forts, however, to establish this phase in real compounds 
largely failed. The most promising candidates so far are 
the hexaborides,^ where the unexpected ferromagnetism 
may be interpreted in terms of a doped excitonic insu- 
lator,'^'^ and TmSeo.45Teo.55, where the unusual thermal 
diffusivity may be due to a superfluid exciton state.^ 

We focus on TmSeo.45Teo.55- In a series of experiments, 
Wachter and coworkers compiled strong evidence for ex- 
citonic phases in the vicinity of the pressure-induced SC- 
SM transition in TmSeo.45Teo.55-^"^ In particular, the 
anomalous increase of the electrical resistivity in a nar- 
row pressure range around 8 kbar indicated the appear- 
ance of a new phase below 250 K. Wachter and coworkers 
suggested that this new phase might be an "excitonic in- 
sulator", and, assuming the energy gap Eg to decrease 
with increasing pressure, constructed a phase diagram for 
TmSco.45Teo.55 in the Eg ~ T plane. ^ Later they found in 
the same material a linear increase of the thermal diffu- 
sivity below 20 K and related this to a superfluid exciton 
state. ^ Both excitonic phases are located on the SC side 
of the SC-SM transition in TmSeo.45Teo.55. 

In our previous work^ we presented a theoretical 
analysis of these astonishing experimental findings. We 
calculated for a model mimicking the situation in 
TmSeo.45Teo.55 the phase boundary Tc{Eg) for an ex- 
citonic insulator and investigated the composition of its 
environment. Our main conclusions were: (i) The phase 
boundary constructed from the resistivity data does not 
embrace the excitonic insulator. Instead, it most proba- 
bly traces the exciton-rich region above an excitonic insu- 
lator, where excitons dominate the total density, provide 
abundant scattering centres for the charge carriers, and 
induce the observed resistivity anomaly, (ii) The linear 



increase of the thermal diffusivity below 20 K, on the 
other hand, could signal exciton condensation. The large 
electron-hole mass asymmetry places the excitonic insu- 
lator entirely on the SC side, where it is supported by 
strongly bound, bosonic excitons. A transition to a su- 
perfluid state, similar to the transition in liquid He 4, 
is thus conceivable and could give rise to the observed 
atypical thermal diffusion. 

We based our calculation on an effective-mass, 
Wannier-type model for electrons in the lowest conduc- 
tion band (CB) and holes in the highest valence band 
(VB), with an indirect energy gap separating the two 
bands, a multiplicity factor to account for three equiva- 
lent CB minima arising from the fee crystal structure of 
TmSeo.45Teo.55, and a statically screened Coulomb po- 
tential operating between the constituents of the model 
(in contrast to the quasi-static screening model used be- 
low). Because of the limitations of the model, which, 
for instance, neglects lattice effects, we could not obtain 
complete quantitative agreement between theory and ex- 
periment. Yet, our results indicate that Wachter and 
coworkers^"^ may have seen an exciton condensate. 

Clearly, within an effective electron-hole model quanti- 
tative agreement between experiment and theory cannot 
be expected. Instead of augmenting the model with ad- 
ditional degrees of freedom, we still want to stay in the 
present paper within its bounds and discuss the possi- 
bility of an electron-hole liquid^ which competes with 
excitonic phases in materials where the band structure 
is multi-valleyed. The results will further benchmark the 
model and thus guide its modification. 

2. Method 

The thermodynamics of an electron-hole liquid is 
contained in the electron-hole pair chemical potential 
^eh {fs , T) which depends on the electron- hole pair den- 
sity parameter Ts and the temperature T.^ In our case, 
^eh = —Eg, with Eg the energy gap fixed by the ex- 
ternal pressure. Therefore, we first have to determine 
rs{Eg, T) and then convert this relation into —Eg(rs,T). 
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Fig. 1. The band structure of the effective-mass, Wannier-type 
model. We measure momenta k for the CB electrons and VB 
holes from the respective extrema of the bands which, in reality, 

are separated by w = K/2 where K is the reeiprocal lattice vec- 
tor connecting the F with the X points of the Brillouin zone for 
a fee lattice. The energy gap Eg varies continuously through zero 
under pressure. Solid (dashed) parabola denote the bare (renor- 
malised) energy bands. Note, the auxiliary clicinical potentials 
/ij are defined with respect to the shifted band edges. The factors 
Qi account for spin and the number of equivalent valleys. 



Comparing —Eg{rs,T) with the exciton binding energy 
and analysing for fixed temperature its rg-dependence 
enables us to determine the region in the Eg — T plane 
where an electron-hole liquid may be favoured over exci- 
tonic phases. We calculate —Eg{rs, T) in the quasi-static 
plasmon-pole approximation.^*^ 

The schematic band structure used in Ref.* to model 
TmSeo.45Tco.55 is shown in Fig. 1. Within the quasi- 
static plasmon-pole approximation the CB electron [i = 
1) and VB hole (i = 2) self energies, leading, at finite fill- 
ing, to the renormalization of the energy bands (dashed 
parabola in Fig. 1), can be written as 

E,(k) = Sf (k) + Sf(k) , 



where the first term, 

sr(k) = -^y,(k- 



k')nF(ei(k')-/Xi) 



denotes the exchange energy and the second, 
5^f(k) = jE[^«(k')-^o(k')] , 



(1) 



(2) 



(3) 



the Coulomb-hole energy. Here, Vo(q) = 47re^/eo(7^ is the 
Coulomb potential reduced by the background dielectric 
constant £0, K(q) is the statically screened Coulomb po- 
tential to be defined below, ei(k) = k^/2mi [mi are the 
effective masses) , and np (e) is the Fermi function. In con- 
trast to the static approximation used in Ref.,^ which 
replaces directly in the Hamiltonian Vo(q) by 14 (q), the 
selfenergies contain now not only the exchange correc- 
tions but also the Coulomb-hole. 

The screened Coulomb potential in Eqs. (2) and (3) 
is the static limit of the dynamically screened Coulomb 



potential, 

F,(q) = V;(q,0) = yo(q)/e(q,0) 
where the dielectric function is given by 



= 1 



with an effective plasmon dispersion. 



a;(q)^=a;^,[l+(f)Vf?( 
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(4) 



(5) 



(6) 



'Qa' " ^4mi 4m2 

constructed in such a way that sum rules are satisfied. 
When the plasmon-pole approximation is constructed 
from the Lindhard dielectric function, the coefficient 
C = 1. It can be however also used as a fit parameter 
to correct for the overestimation of screening within the 
quasi-static approximation. In that case 1 < C < 4.^° 
The screening wave number. 
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(7) 



(8) 



entering Eq. (6) depend on the CB electron and VB hole 
densities given, respectively, by 



ni=gi'^ «F(ei(k) - m) , 

k 

ri2 = 52EriF(e2(k) - /U2) , 



(9a) 



(9b) 



where the multiplicity factors gi contain the spin and the 
number of equivalent valleys. 

The main effect of the selfenergies Ej is a rigid, k- 
independent shift of the single particle dispersions lead- 
ing to a renormalised band gap Eg. This situation has 
been anticipated in Eqs. (2), (9a), and (9b), where we in- 
troduced auxiliary chemical potentials fii , measured from 
the shifted band edges (see Fig. 1). By construction. 



9 ■ 



(10) 



The bare energy gap Eg, that is, the parameter which is 
experimentally controlled via pressure, is given by 



Eg = Eg 



I]l(0)-I]2(0) 



(11) 



Supplementing Eqs. (10) and (11) with the condition of 
charge neutrality, 

m = n2 , (12) 

we obtain a closed set of equations for the three unknown 
parameters fii and 112, and Eg. 

To construct from Eqs. (10) - (12) the electron-hole 
pair chemical potential —Eg{ra,T), we proceed as fol- 
lows: First, we calculate from Eqs. (12), (10), (9b) and 
(9a), the auxiliary chemical potentials as a function of 
Eg and T. Inserting ^ii{Eg, T) in Eqs. (7) and (8), we ob- 
tain qs[Eg,T) and u>pi{Eg, T), which enables us to deter- 
mine Vs and thus the selfenergies Sj in the whole Eg—T 
plane. Via Eq. (11) we then obtain the important relation 
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Fig. 2. Shifted energy gap Eg as a function of the pressure-tuned 
energy gap Eg and temperature T for TmSeo.45Teo.55 parame- 
ters: a = 0.015, gi = 6, and 92 = 2. The dotted line denotes the 
energy where the excitonic instability would occur at T = and 
the dashed line gives Eg{Eg,T = 0.152) in the static approxi- 
mation. 
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Fig. 3. Electron-hole pair chemical potential —Eg{rs,T) as a 
function of electron-hole pair density parameter rg and temper- 
ature T for TmSeo.45Teo.55 parameters: a = 0.015, gi = 6, and 
32 = 2. The dotted line indicates where the excitonic instability 
would occur at T = and the dashed line gives —Eg{rs,T = 
0.152) in the static approximation. Note, for T < 0.305, the 
isotherms go through a minimum. 



Eg{Eg,T). Finally, we use either Eq. (9a) or Eq. (9b) to 
map, for a given temperature T, the renormalised energy 
gap Eg on to the electron-hole pair density n = ni = n2. 
Using the density parameter — (3/47m • a^j-)^/^, with 
ax the exciton Bohr radius defined in the next section, 
instead of the pair density, we find —Eg{rs, T). 

3. Results 

In the following, we measure energies and lengths in 
units of the exciton Rydberg Rx = l/2ma^ and the ex- 
citon Bohr radius ax = eo/we^, respectively {h = 1). 
The model shown in Fig. 1 is then completely speci- 
fied by the multiplicity factors gi and the effective mass 
ratio a = m\jmi. For TmSeo.45Teo.55 the correspond- 
ing values are g\ = 6, 52 = 2, and a « 0.015.^'^ Due 
to the lack of data about band gap renormalization in 
TmSeo.45Teo.55, we cannot use C in Eq. (6) as a fit pa- 
rameter. Thus, C = 1, which almost certainly overesti- 
mates screening. 

The shifted energy gap Eg as a function of the energy 
gap Eg is shown in Fig. 2 for several temperatures. First 
we note the drastic difference to the static approximation 
(Vb Vs in the Hamiltonian) . Whereas in the static ap- 
proximation Eg decreases monotonously with Eg for all 
temperatures, the quasi-static approximation gives, for 
some E'g-range, a triple-valued Eg for T < 0.305. Be- 
cause of the one-to-one correspondence between Eg and 
the pair density n, this means the system cannot choose 
a unique pair density. As a consequence, there may be 
parts in the -Eg— T-plane where an electron-hole gas coex- 
ists with an electron-hole liquid. Note, the renormalised 
band gap Eg resembles the selfenergy parameter S intro- 
duced before to set up a Hartree-Fock description of the 
gas/liquid transition in a generic electron-hole plasma. 

The corresponding pair chemical potential —Egijs^ T) 
is shown in Fig. 3. Various subtle points should be no- 
ticed before we can attempt an interpretation of the data. 
First, —Eg{rs, T) does not asymptotically merge with the 
exciton Rydberg for 00. Thus, the results in the ex- 



treme dilute limit are flawed, not only in the quasi-static 
approximation (solid lines) but also in the static approx- 
imation (dashed line). This is a well-known shortcoming 
of random-phase-approximation based screening models, 
which neglect multiple-scattering between charged car- 
riers. Including multiple-scattering (ladder diagrams) in 
the quasi-static (static) approximation would force the 
chemical potential —Eg(rs,T) to approach for Vs 00 
the exciton Rydberg from below (above). Second, the 
densities of interest are rather low [vs > 1). Corrections 
due to multiple scattering are therefore not necessarily 
small and should be included, a formidable task beyond 
the scope of the present paper. 

With these caveats in mind we now extract the phys- 
ical implications of the data shown in Fig. 3. As ex- 
pected from the E'g-isotherms, the model shown in Fig. 1 
yields for TmSeo.45Teo.55 parameters a chemical poten- 
tial which shows the typical van-der-Waals behaviour: 
Above a critical temperature {Tcrit ~ 0.254), the chem- 
ical potential increases monotonously with density (re- 
call: Vs ^ 1/n^/^). Below the critical temperature, the 
chemical potential has a local maximum and a local min- 
imum. This can be most clearly seen from the T — 0.203 
isotherm. A Maxwell construction would thus give the 
equilibrium chemical potential and the coexistence re- 
gion between the gaseous and the liquid phases. An- 
ticipating however the correct dilute limit, a Maxwell 
construction cannot be made and the minimum sim- 
ply indicates a jump in the pair densities, that is, for 
large enough energy gaps the pair density n drops to 
zero.^^ Only when the exciton binding energy is below 
the minimum of the chemical potential, a gas/liquid tran- 
sition and with it a coexistence region could appear. 
In that case, excitons would be also stable in a certain 
range. Because of the experimental evidence for excitonic 
phases, we expect this situation to be actually realised 
in TmSeo.45Teo.55- 

The minimum of the pair chemical potential 
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Fig. 4. Phase boundciries obtained for TmSeo.4BTeo.65 parame- 
ters: a = 0.015, PI = 6, and p2 = 2. The thick dotted line is the 
50%-contour of the bound state fraction 7.® Beyond the dashed 
line excitons disappear because of the Mott effect and the grey 
region denotes the excitonic insulator (El).** At the full diamond 
on the horizontal thin dotted line the linear increase of the ther- 
mal diffusivity has been experimentally found'' and at the full 
circle n''^ 1.3 X 10^^ cm~^* in good agreement with the exper- 
imentally estimated exciton density 3.9 x 10^^ cm~^.^ The open 
circles are the points where —Eg{rs,T) goes through the mini- 
mum. The energy gap of TmSeo.4BTeo.5B at ambient pressure is 
given by the vertical thin dotted line. 

—Eg{rs,T) is below the exciton Rydberg. Therefore, tak- 
ing the correct ^ 00 Hmit of —Eg{rs,T) into ac- 
count, our results for the pair chemical potential suggest 
that in the Wannier-type model excitons and thus exci- 
tonic phases (exciton gas, excitonic insulator) are unsta- 
ble against a metallic phase. Anticipating, in addition, 
the decrease of the exciton binding energy with increas- 
ing density, the metallic phase would be in fact favoured 
over most of the relevant parameter range discussed in 
Ref.^ and in stark contrast to the experimentally ob- 
served increase of the resistivity. 

This can be seen in Fig. 4 which combines the points 
where —Eg (r^ , T) is minimal with (i) the phase boundary 
of the excitonic insulator, (ii) the 50%-contour of the 
bound state fraction 7 = n\/{n{ + nj) {n( and n\ are 
the free and the bound part of the total density obtained 
from a T- matrix calculation of the total density) , and (iii) 
the Mott line beyond which excitons disappear because 
Pauli blocking and screening lead to a vanishing exciton 
binding energy. All of it has been calculated in Ref.® 
within the static screening model. We expect however 
the quasi-static approximation to give similar results. 

One of our main conclusions in Ref.^ was that the re- 
sistivity anomaly, that is, the experimentally observed 
increase of the resistivity in a narrow pressure range, is 
most probably due to c;xciton-elcctron and cxciton-hole 
scattering which would be very efficient above an exci- 
tonic insulator. Indeed, the temperature range, where the 
50%-contour is located, coincides roughly with the tem- 
perature range where the anomaly has been detected. 
In contrast, the analysis of the pair chemical potential 
performed in this paper suggests excitons in that region 
to be unstable against a metallic electron-hole plasma. 
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An increase of the resistivity is however hard to recon- 
cile with a metallic phase. We expect therefore excitons 
in TmSeo.45Teo.55 to be stabilised by band structure ef- 
fects beyond the model shown in Fig. 1 and additional 
scattering processes. That the effective-mass, Wannier- 
type model may not be quite adequate can be also seen 
from the fact that the experimentally determined exci- 
ton Rydberg for TmSeo.45Teo.55 is » 75 meV whereas 
the energy gap at ambient pressure is ~ 135 meV and 
thus only w 1.8 • Rx (vertical thin dotted line in Fig. 4). 

In analogy to the hexaborides,^^ we expect coherent 
inter- valley scattering of excitons to increase the exciton 
binding energy. Additionally, exciton-phonon scattering 
should also increase the binding energy of an exciton in 
TmSeo.45Teo.55 because the narrow VB makes the exci- 
ton very susceptible to phonon dressing. In fact, there 
is plenty of experimental evidence for strong exciton- 
phonon coupling in TmSeo.45Teo.55.^ Both scattering 
processes should be studied in detail, ideally in a model 
which avoids the effective mass approximation and keeps 
the full band structure. If it turns out that in a more re- 
alistic model the exciton binding energy is still above the 
chemical potential of the metallic phase, excitonic phases 
must be ruled out in TmSeo.45Teo.55. On the other hand, 
if the exciton binding energy is pushed below the mini- 
mum of the metallic chemical potential, and that is not 
unrealistic considering the energy gain coherent inter- 
valley scattering alone can yield if the band structure 
permits it,^^ excitons are stable, at least in some pres- 
sure range, and the excitonic instability leading to an 
excitonic insulator could take place in TmSeo.45Teo.55 
before an electron-hole gas-liquid transition interferes. 
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